Int. J. Adv. Sci. Eng. Vol.7 No.3 1849-1853 (2021) 1849 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


Numerical Simulation of Novel Hybrid Solid 
Desiccant Bed for Drying 


Selvaraji Muthu}, N. Sekarapandian2* 
1 Department of Thermal and Energy, School of Mechanical Engineering, Vellore Institute of Technology, 
Vellore Campus, Tiruvalam Road, Katpadi, Vellore, Tamil Nadu 632014, India 


2Department of Thermal and Energy, School of Mechanical Engineering, Vellore Institute of Technology, 
Vellore Campus, Tiruvalam Road, Katpadi, Vellore, Tamil Nadu 632014, India 


ABSTRACT: In this paper, a novel design of hybrid stationary packed solid desiccant bed is proposed and 
numerical simulation of the performance for drying of gases is presented. The pro-posed hybrid packed beds is 
consists of a mixture of solid desiccants of different types having both fixed and varying particle diameter 
distribution along the vertical axis direction. Generally in the beds of this type, the drying time taken for the bed 
saturation is hugely related to the input condition of the gases. Also, the adsorption isotherm curves of different 
desiccant types have a strong effect on the exit conditions of the gases dehumidified/dried using these beds. In 
order to simulate the transport phenomena associated with the hybrid packed solid desiccant bed application, an 
in-house CFD code using finite volume method is developed. This CFD code has the capability to model the 
simultaneous heat conduction and moisture diffusion taking place in the bed during the gas transport. The 
numerical results of the base model design are validated against the experimental data available in the literature. 
The estimated exit temperature and moisture content of the gas are in good agreement with the literature data. 
Further, the validated model is used to simulate for several other novel packed bed designs, the coupled heat and 
mass transfer in adsorption and desorption phases under wide operating conditions. The results of interest 
includes the assessment of improvement in performance of dehumidifying/drying of gases in terms of pressure 
dew point, relative humidity at exit and the prediction of the desorption phenomena for the regeneration of the 
beds for next cycle of operation. Finally, the paper summarizes the benefits of this novel hybrid design and 
highlights the potential design cost benefit. 
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INTRODUCTION improve the adsorption and desorption phenomena 


The stationary hybrid fixed desiccant bed is filled with 
solid desiccant particles as shown in Figure 1. The 
moist air is supplied to the bed and dry hot air exits at 
the outlet during the adsorption process [1]. During 
the adsorption process water-vapor molecules are 
transferred from the bulk air-flow to the desiccant 
surface due to vapour pressure difference, and the 
direction of mass transfer is inversed during the 
desorption process and the desiccant bed gets 
regenerated. As the moisture vapour pressure on the 
surface of the desiccant particle is highly influenced 
by the desiccant temperature and the moisture 
content on the desiccant, the heat of adsorption is 
dependent on the mass transfer rate, which is mainly 
depended by the moisture vapour pressure on the 
desiccant surface. Hence, this issue of adsorption and 
desorption in the hybrid fixed desiccant bed must be 
treated as a coupled transient heat and mass transfer 
phenomena, which are numerically simulated [2] and 
experimental measure of the performance of the beds 
[3 & 4]. This adsorption and desorption phenomena is 
analyzed for rotary desiccant wheel application as 
well [5-7]. 


A novel hybrid packed beds [8] are proposed to 


during the drying application. The hybrid desiccant 
beds [9-12] show the capability to achieve a further 
low dew point temperature for the drying 
applications. 


The concept of novel hybrid desiccant packed bed 
proposed along with conventional Silica gel (SG) and 
Molecular Sieves (MS) beds for performance 
comparison purpose is shown in Figure 1. 


° SG desiccant bed 
° MS desiccant bed 
° Hybrid bed of 1/2 SG and 1/2 MS desiccant 


In this hybrid desiccant bed, the half of the bed 
volume at the bottom is filled with SG desiccant and 
the remaining half of the portion at the top of the 
tower is filled with MS desiccant (Zeolites). By this 
novel concept, the higher adsorption capacity of the 
silicagel desiccant is utilized to adsorb the maxi-mum 
moisture at the entry portion of the tower and then 
the remaining moisture content in the air is effectively 
adsorbed using molecular sieves desiccant, as 
molecular desiccant is having much higher adsorption 
capacity than silicagel at lower relative humidity 
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conditions. This is a unique advantage of this hybrid 
packed bed as compared to the conventional beds of 
entirely using either silicagel or molecular sieves. 


The objective of this paper is to present the novel 
design of hybrid stationary packed solid desiccant bed 
to achieve superior DPT while retaining the same size 
and volume for the conventional SG and MS bed for 
comparison purpose. Also the CFD code is to simulate 
the performance of the packed bed and validating the 
results with experimental results from literature and 
then to present to the potential benefit of using hybrid 
bed in the place of conventional packed beds. 


Bed with Silica gel’ Bed with Molecular 
Activated Alumina 


Hybrid Bed 
Sieves #1 


Fig. 1 Hybrid stationary packed solid desiccant bed 


Governing equations 

The governing equations in the form of partial 
differential equations are derived [3, 4] based on the 
conservation of mass/moisture and conservation of 
energy/heat in the air and in the desiccant. 


The mass/moisture conservation of the air in the solid 
desiccant packed bed is given by partial differential 
equation as below: 


OW ka 
= m (a, @,) 
ot (-e)p, (1) 


Mass/moisture conservation in desiccant for the solid 
desiccant packed bed is given by partial differential 
equation as below: 

00, G, 0a, 2 k,,a aa (2) 

Ot gp, OY €&, 

Heat energy conservation in the solid desiccant 
packed bed air is given by partial differential equation 
as below: 


20, G, 00, k,a 
= + 


a 


(6, 8, ) 


at p,é by p,C,é (3) 
+Gknt (a, 0,)(0, -,) 
C,Paé 
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Heat energy conservation in the desiccant for the solid 
desiccant packed bed is given by partial differential 
equation as below: 


08, Ag C8 ® ka 
| a : | m (a, Q,) 
ot (l-e)p,c, Oy oc, d-e)p, 
+H“ _6,-0,) 
(l-€)p,c, (4) 


The initial, periodic and boundary conditions are 
referred from literature [2] to compare the numerical 
results [2,3]. 


Numerical modeling and methodology of solving 
governing equations 


The governing partial differential equations are 
discretized by using finite volume method under 
implicit scheme to get set of linear algebraic equations 
and these linear equations are solved with the 
boundary conditions and input conditions using the 
CFD code developed in FORTRAN 95, as shown by the 
flow diagram in Figure 2. 


Solving Equation-1 


Solving Equation-2 


Soiree Eaton’ 
Calculation of properties 


Meshing Control volumes 


Fig. 2 Modeling and methodology of solving 
governing equations 


Validation 

The simulated results of present model are compared 
with numerical results [1 & 2] and experimental [3] 
results from the literature. The parameters of packed 
bed are shown in Table-1. 


The numerical model shows a finer agreement with 
experimental data for the exit air temperature and 
moisture content of process air during adsorption, as 
shown in Figure 3. 
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RESULTS AND DISCUSSION 


The moisture content and air temperature for the bed 
sections across the length for various inlet moisture 
content from 5-20 g/kg are predicted by numerical 
model, as shown in Figure 4. It is evident that the 
process air temperature is directly depending on the 
moisture content removal, which is in turn related to 
the inlet moisture content of process air. 


Table 1. Parameters of the packed beds 


Diameter of the 130 

packed desiccant a 

beds 

Diameter of the 388mm 

desiccant particles 

Height of the 77.55mm 

desiccant packed 

Desi 4 Molecular 

Belecani ype tse Sieves(MS) , Silicagel 

(SG) 

Temperature of 23.3 0C 

Process air inlet : 

Moisture content of 18 g/kg of dry air 

Process air inlet 

Relative humidity of 100% 

Process air inlet 

Porosity of the 0.4 

desiccants 

Ratio of air flow é 

passage area to bulk 81.6% 

area 

Maximum 0.35 for SG, 0.25 for 

adsorption capacity, MS 

Operating time of 

the packed bed, t 1800 sec 
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os 


° Ta Pesaran rant Experimental 
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wke 


4 Ya Pesaran run! Expermnental 


Ya Ramzy rua l-Numerical 


}, ~ 388 mm L = 77.5 oom, v O21 ms 


Exit temperature of process air ( Ta ), deg C 


Moisture cotent of process air at exit ( Ya), 


W, = 41.7Z, Ty 23.3'C, £1800 se 
ke 


0.00 0.20 0.40 0.60 0.80 1.00 


Relative time of adsorption cycle 


Fig. 3 Validation of numerical results with 
experimental results 
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The corresponding moisture removal capacity (MRC) 
and relative moisture removal efficiency (RMRE) are 
plotted for the bed sections across the length for 
various inlet moisture content from 5- 20 g/kg is 
shown in Figure 5. It is seen that MRC is directly 
depending on the inlet moisture content of process 
air, however, the RMRE is inversely proposal to inlet 
moisture content. 


OTg@Yg=20 g/kg 
=Tg@Yg=5 g/kg 
© Yg=10 g/kg 


XTg@Yg=15 g/kg 
* Yg=20 g/kg 
+Yg=5 g/kg 


OTg @ Yg=10 g/kg 
AYg=15 g/kg 


oOo wow fF mn WB Nv 
mM mM mM MN Nn MN 


Exhaust Temperature of air, C 
n 


Moisture content of air at exhaust, g/kg 


Nn 


1 
mn 


0.0 0.2 0.4 0.6 0.8 1.0 


Desiccant Fixed bed operating time 


Fig. 4.Calculation of exit temperature for different 
inlet moisture and temperature 
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Fig. 5.Calculation of RMRE and MRC for different 
inlet moisture and temperature 


Figure 6 shows the relative humidity in % for the bed 
sections across the length. It is seen that the moisture 
removal is very strong with the shorter length of bed, 
and the instantaneous relative humidity drop is very 
significant (that is from 100% at inlet, it is dried to 
5.5% within 50% of the bed length).So in order to 
display the further drop of moisture contents, a 
logarithmic scale is used for Y-axis in all the RH 
graphs, a lowest value of 0.4% is able to be displayed 
clearly. 
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For the plots of DPT values achieved over the bed 
length, a conventional normal scale is used, as it can 
be seen in Figure 7. At the outlet a lowest DPT of 
about - 39.6°C using the hybrid bed #1 (50% : 50% = 
SG :MS) is shown in Figure 8 and Table 2, at the same 
time, the conventional MS bed gives a DPT of -39.01°C. 


=) 
S 
i) 
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=SG:MS 
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=SG:MS 
amie (1)%:100 % 
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= 
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Fig. 6.Hybrid packed desiccant bed 
performance (RH) for 18 g/kg 
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Fig. 7 Hybrid packed desiccant bed performance 
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Fig 8 Performance of hybrid packed bed 
at 18 g/kg inlet 
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CONCLUSIONS 


The numerical model shows a finer agreement 
with experimental data. The effect on output moisture 
content and temperature of air are predicted under 
the various input moisture contents. 

The effect on output moisture content and DPT of 
air are predicted under different input moisture 
contents for the conventional SG, MS beds and hybrid 
bed. 

The Hybrid bed with 50% of bed length with SG 
desiccant and remaining 50% of bed length with MS 
desiccant is superior to the conventional MS bed, by 
which the cost of hybrid bed can be lower than that of 
MS bed, without compromising the performance. 


Table 2 DPT at outlet of the hybrid desiccant bed 


Highest achieved Dew Point Temperature 
(°C) 
Inlet 
SGbed | Hybridbea | MS 
; bed 
Moisture 
18 g/kg -17.17 -39.60 -39.01 
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Definitions/Abbreviations 


Db Diameter of the bed (m) 

Dp Diameter of the desiccant particle (m) 

Cp specific heat at constant pressure (kJ/kg K) 
Dx Knudsen diffusivity (m2/s) 

Do ordinary or molecular diffusivity (m?2/s) 

Ds surface diffusivity (m2/s) 

k thermal conductivity (W/m K) 

Ky Convective Mass Transfer coefficient (kg/m? 
h Te heat transfer coefficient 

(kW/m?2k) 

L Length of the bed (m) 

Le Lewis number 

Nu Nusselt number 

P pressure (Pa) 

Pa Atmospheric pressure (Pa) 

st heat of adsorption J/kg desiccant 
r pore radius of the desiccant (m) 

RH Relative humidity (%) 

Sh Sherwood number 

- temperature (K) 
t Time 

u velocity (m/s) 

WwW adsorption in kg water vapour/kg desiccant 
bg Humidity ratio, kg water vapour/kg dry air 
Z axial direction 

Greek 
symbols 
p density (kg/m3) 
€ porosity 
Subscripts 

0 initial state 

a air 

d desiccant 

p process air 
r Regeneration air 
s Solid desiccant 
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